Topology and existence of 3D anisotropic filamentary 

kinematic dynamos 

by 

L.C. Garcia de Andrade 

s 

o 

Departamento de Ffsica Teorica - IF - Universidade do Estado do Rio de Janeiro-UERJ 

Rua Sao Francisco Xavier, 524 
^ ! Cep 20550-003, Maracana, Rio de Janeiro, RJ, Brasil 

Electronic mail address: garcia@dft.if.uerj.br 



> 

O 

00 



Abstract 



Curvature and helicity topological bounds for the magnetic energy of the streamlines mag- 



CD 

netic structures of a kinematic dynamo flow are computed. The existence of the filament 
dynamos are determined by solving the magnetohydrodynamic equations for 3D flows and the 

Q-f 

solution is used to determine these bounds. It is shown that in the limit of zero resistivity 

& ■ 

filamentary dynamos always exists in the isotropic case, however when one takes into account 
that the Frenet frame does not depend only of the filament length parameter s, (anisotropic 
^ ■ case) the existence of the filamentary dynamo structure depends on the curvature in the case 

of screwed dynamos. Frenet curvature is associated with forld and torsion to twist which allows 
us to have a sretch, twist, and fold method to build fast filament dynamos. Arnold theorem 
for the helicity bounds of energy of a divergence- free vector field is satisfied for these stream- 
lines and the constant which depends on the size of the compact domain MCR 2 , where the 
vector field is defined is determined in terms of the dimensions of the constant cross-section 
filament. It is shown that when the Arnold theorem is violated by the filament amplification of 
the magnetic field structure appears, the magnetic field decays in space. PACS numbers: 

02.40.Hw:Riemannian geometries 
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I Introduction 



The topology and geometry of hydrodynamical and MHD dynamos have been mainly developed 
by Arnold and Khesin [1] and by Childress and Gilbert [2] using non chaotic flows and the twist, 
stretch and fold technique developed by Moffatt [3] to investigate fast dynamos. Dynamos 
in chaotic flows have been developed lately by Thiffault and Boozer [4]. On the other hand 
twisted filamentary magnetic structures have been important in plasma and solar physics [5] in 
the investigation of electric carrying-current loops. Though we know [2] that planar dynamos 
do not exist they amplify the magnetic fields. It is also important to recognize when we have a 
dynamo or anti-dynamo from the magnetic profile. Since we only have planar dynamos in very 
specific circunstances, when we have folds of filaments for example, we here consider three 
dimensional 3D dynamos where the twist can be associated with torsion of the filament and 
curvature to fold processes. In this paper we consider the topology and geometrical bound 
obtained from the solution of magnetic MHD equations with a dissipative term in the absence 
of electric potential. The scalar MHD equations to be solved are obtained by expressing the 
MHD equations in the Frenet frame anisotropic basis where the basis depend upon not only 
the position of the frame along the filament but also of the normal and binormal direction 
of the frame as well as the time. The filament are considered to be planar with constant 
curvature. The existence of anti-dynamos where the magnetic field is amplifies, as in the case 
of galactic magnetic fields, depends upon the relation between the right (on the left) handness 
of curvature of the filament and the normal direction to the filament which is assumed to 
move along the binormal direction such as some vortex filaments. Actually this is not the 
first time Frenet curvature is used to investigate dynamos, in 2001 Schekochihin et al [6]used 
the statistical Frenet curvature in the small-scale magnetic fields in kinematic dynamos. This 
paper is organized as follows: In section 2 we compute the filamentary solution of the magnetic 
filament. In section 3 we compute the energy bounds from helicity and curvature and examine 
the relation with the existence of the filamentary dynamos. The conclusions are presented in 
section 4. 
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II MHD scalar equations for kinematical dynamos 



Let us now start by considering the MHD field equations 

V.S = (II.l) 

d 

—B-X7x[uxB]-eX7 2 B = (II.2) 
where u is a solenoidal field while e is the diffusion coefficient. Equation (11.2) represents the 

— * 

induction equation. The magnetic field B is chosen to lie along the filament and is defined 
by the expression B = B(s,n)t and u = ub, chosen by the similarity with vortex filaments, 
is the speed of the flow. The remaining coordinate n is orthogonal to the filament all along 
its extension, and the arc length s measures distances along the the filament itself. The 
vectors t and n along with binormal vector b together form the Frenet frame which obeys the 
Frenet-Serret equations 

? = /era (II.3) 

n' = -nt + rb (II.4) 

b' = -rn (II.5) 

the dash represents the ordinary derivation with respect to coordinate s, and n(s,t) is the 
curvature of the curve where k = Rr 1 . Here r represents the Frenet torsion. We follow the 
assumption that the Frenet frame [7] may depend on other degrees of freedom such as that 
the gradient operator becomes 

->d d -*d 

v = t & + % + ^ (II6) 

The other equations for the other legs of the Frenet frame are 

d 

—t = 9 ns n+[n b + r]b (11.7) 



dn 

d_ 

dn 



n = -6 ns t- [divb)b (11.8) 



& - 

—b = -\Q b + r]t- (divb)n (11.9) 
on 

—t = 9 bs b-[a n + r]n (11.10) 



d _ 

—n=[tt n + T]t-K+(divn)b (H.ll) 

d - 

—b = -9 bs t - [k + (divn)]n (11.12) 

Another set of equations which we shall need here is the time derivative of the Frenet frame 
given by 

t = [n'b-KTn] (11.13) 

k=KTt (11.14) 

b=-Kt (11.15) 

A long and straithforward computation .specially due to the computation of V 2 A, where the 

vector potential A = A(t,s,n) in principle. Substituting these equations for the dynamics of 
the Frenet frame leads to the scalar MHD expressions 

d t A = -d s <P + [d\A - A(6 r J - k 2 )] (11.16) 

- KT A = -uB + e[2d n A + (Q s + T)6 ns A] (11.17) 

-9 bs A = e[2d n AVL s + tt 2 A] (II. 18) 

where « is the Frenet curvature of the streamlines. These equations have already been 
simplified by using the relations 

VxA = B (11.19) 
which yields the following differential scalar equations 

B = -A[n b + r] (11.20) 

d n A + KA = (11.21) 

A(fi n + r) = (11.22) 

Where the f2's represent the abnormalities of the streamlines of the flow. Here 

0ns = n-^-t (11.23) 

on 
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When the £l s vanishes we note the geodesic streamlines are obtained. As we shall see below 
here we are not consider geodesic flows dynamos. By considering planar flows where torsion 
vanishes and the gauge condition 

V.I+J^ = (11.24) 

This equation can be expressed as 

d a A+[6 na + 9 ba ]A = (11.25) 
Now by considering that A does not depend on the coordinate s this expression reduces to 

[Ons + e bs ]A = (11.26) 

which reduces to 6 ns = —9 bs . By making use of this expression and the assumption that = 
one simplifies the MHD scalar equations to 

d t A = [d\A - A(9 n 2 - k 2 )] (11.27) 

-k tqA = -uA[n b + r ] + e[2K + (Sl a + T )9 ns }A (11.28) 

Where we have used the hypothesis of helical or screwed dynamos where torsion and curvature 
coincides and are constants and equal to r and k . Simple algebraic manipulation of these 
equations, in the limit of zero resistivity ,or e = reduce them to 

k 2 + uk + Vt b u = (11.29) 

which is a second order algebraic equation, and 

d t A = [2n 2 - 9 ns 2 ]A (11.30) 

d n A = n A (11.31) 

The discriminant of the algebraic equation is assumed to vanish or A = u 2 — AVL b u = 
which yields the solution u = — AVt b . A simple solution of the system can now be obtained 
by separation of variables A(t,n) = H(n)T(t), substitution of this product in the equations 
yields 

A = A e lxt ~ Kon] (11.32) 



which by the relation between the fields B and A yields 

B = A (n b + T )^ xt - K ^ (11.33) 

where A := [2k 2 — 9 ns 2 ]. Both results shows that the filament magnetic field can only be 
maintained in time when A > or 2k 2 > 6 ns 2 otherwise the magnetic field decay very fast in 
time and the filament is actually an antidynamo. Note that when the filament is isotropic in 
the sense the frame only depends upon the arc length s then the filaments are always dynamos 
since then 0„ s = 0. 



Ill Topology bounds the energy from magnetic he- 
licity 

The magnetic energy of a divergence vector field on a compact planar domain M reads 

E B = — I B 2 dV (111.34) 
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substitution of the formulas obtained in the last section yields 

2 

Eb = jJ A 2 [Q b + r 2 ]ds (111.35) 

where we consider that the magnetic filaments leads constant cross-section area as S = na 2 . 
In terms of the vector potential component A n expression (111.35) reduces to 

T a 2 

Eb = ^[B 2 ] (111.36) 

where L — J ds is the length of the filament. Now the Arnold theorem states that a divergence- 
free vector field a defined a compact manifold M , the helicity have an upper limit exactly 
given by the modulus of the magnetic helity 



H = J a.VxadV (111.37) 
where the topological bound is given by [?] 

E B > C\H\ (111.38) 



where C is a positive constant which depends upon the size and form of the compact manifold. 
Let us now apply the Arnold's theorem to our solution and analyze the implications for the 
anti-dynamo problem. First of all we must compute the magnetic helicity in the example given 
here. This yields 

\H\ = na 2 L[n b + T ]A 2 e [xt - K ° n] (111.39) 
Thus to examine the comparison with the Arnold's theorem we have 

C = ^ = ^ b + T ] (HI.40) 

which is the Arnold's constant C which clearly fulfill the geometrical requirements of Arnold's 
theorem. One must notice that when the dynamo conditions discussed in the previous section is 
fulfilled the Arnolds theorem is also obeyed, however when dynamos existence is not possible 
Arnold's theorem seems to be violated and an lower bound for the energy is not obtained 
anymore. 

IV Conclusions 

While the existence of planar MHD solutions does not warrant the existence of dynamos (with 
the possible exception of the presence of fold of filaments) it is enough to warrant the existence 
of the amplification of magnetic fields. Unfortunately a simple way, though not the only, of 
getting fold of the filaments is by introducing Frenet torsion in the filaments but this would 
violate the condition of planarity. In conclusion, the investigation of kinematical anti-dynamos 
and dynamos filamentary MHD in generalised Frenet frame shows that is possible to test the 
Arnold's theorem against the dynamo conditions on Frenet curvature and topology comparing 
it with the magnetic helicity of streamlines. We pretend in near future to investigate this 
theorem for a more general class of fluids where the topological numbers of twist and writhe 
can be computed for other more complex filamentary geometries. 
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